■~^23 

N78-32489 


ACCURACY OF RESULTS WITH NASTRAN MODAL SYNTHESIS 

D. N. Her ting 
Universal Analytics, Inc. 


SUMMARY 


A new method for component mode synthesis has been developed for installa- 
tion in NASTRAN Level 17.5. An introduction and summary of the method was 
presented at the 1977 NASTRAN Colloquium [Ref. 1], but actual results were 
unavailable at that time. This paper serves as a continuation to Reference 1 
by presenting results obtained from the new methoa and comparing these re- 
sults with existing modal synthesis methods. 


INTRODUCTION 


The modal synthesis system developed by Universal Analytics, Inc. (UAI) for 
NASTRAN is a new development which provides for the benefits inherent in 
existing methods but eliminates the restrictiJns and computational drawbacks 
associated with other methods. In Reference 1 it was postulated that the new 
method was sufficiently general to duplicate the results of other, more re- 
stricted, methods simply by choosing different types of normal modes or 
vector recovery procedures. The test problems described herein have been 
selected for direct comparison with other published results. The goal of 
this effort was to determine the relative accuracy of the UAI method with its 
different options. 

The use of structural modes as generalized degrees of freedom in dynamic models 
originated in the analog computer field where structures were combined with 
aeroelastic and control system models. The first applications to digital com- 
puters were simple extensions of the analog techniques. This so-called clas- 
sical approach proved both highly restrictive and limited in accuracy. Many 
different approaches have been developed in recent years having increased 
accuracy and more generality in solving large-order structure dynamics pro- 
blems . 

Although the previous methods used in component mode synthesis vary consider- 
ably in both approach and application, they may be grouped into two distinct 
categories. The first category contains all of the methods using a Rayleigh- 
Ritz approach in which the component degrees of freedom represent the deflec- 
tions of normal modes and static deflection shapes. The second category 
contains methods in which the component degrees of freedom are actual 
physical displacements plus a set of modal coordinates. Here, the classical 
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method has been improved by adding flexibility coefficients to the matrices to 
account for the effects of a truncated set of modes. 

The theoretical development of the NASTRAN modal synthesis system is being 
issued in Reference 2, the Level 17.5 Theoretical Manual. In this develop- 
ment the "residual flexibility" approach is used as a starting point but the 
end results become very similar to the Rayleigh-Ritz approach. The new 
method, in effect, is related to both categories of modal synthesis and shows 
that the differences between chem are more related to computational proce- 
dures than in theoretical basis. 

The two test problems described below were selected for comparison with several 
advanced mode synthesis methods. The problems also provide a comparison of 
the various options that will be available In the NASTRAN system which are 
summarized below: 

1. The boundary conditions used to obtain component modes are not restricted. 
Free, constrained, and partially free modes may be used. 

2. Inertia relief displacement shape functions may be Included as degree's of 
freedom as a user option. These provide for exact static response of 
free bodies and more accuracy for low frequency response. 

3. In the vector recovery process, after a system solution has been obtained, 
a "mode acceleration" procedure which calculates "Improved" displace- 
ments is available. 

4. A full set of error check procedures are available to assess the accuracy 
of the results. These Include printout of the equilibrium forces, energy 
checks of truncated modes, and direct evaluation of the participation of 
the modal coordinates . 

The test problems and their results are summarized below, followed by a 
summary of the conclusions which follow from the evaluation of the tests. 


NOMENCLATURE 


u - Physical Displacement 
G - Guyan Reduction Transformation Matrix 
K - Stiffness Matrix 
M - Mass Matrix 
it - Length 

x - Spacial Coordinates 
T - Kinetic Energy 
V - Potential Energy 
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e - Error Ratio 

£ - Generalized Displacement of a Mode 

P - Density 

<|i - Eigenvector 

10 - Radian Frequency 


EXACT ROD PROBLEM 


A convenient test problem Cor modal synthesis evaluation was used by Rubin 
[Ref. 3] to compare various methods, including his own new method. The prob- 
lem, illustrated in Figure 1 consists of a single rod with extenslonal motion. 
Rather than solve the problem with finite elements, a set of closed form in- 
tegral solutions may be used to obtain the modal synthesis matrix coefficients. 
In effect, the results will simulate a problem with an infinite number of 
elements. This procedure will elminate the finite element errors and will 
allow analysis of errors resulting only from the modal synthesis formulation. 

The problem solved by Rubin uses the free-free modes of the rod to fomulate 
a component mode substructure. The solution matrix is then constrained to 
obtain cantilever modes. If the end degree of freedom were included in the 
normal formulation It could be attached to another structure directly. The 
errors in the results will occur because the sine wave solutions for the 
cantilever rod must be approximated by the dissimilar cosine waves of the 
free rod. 


In the UAI method, the displacement shapes are 
Static Displacements 


u^(x) = 
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Normal Modes: 


V x) " ( +tt - '■■♦biA 

- ( cos - l) e k k - 1,2,...“ (3) 

After normalizing Che units, the total displacement at any point on the rod is: 

u(s > • “i + l( x - fr) 5 o + ^(“ s TT - l ) 5 k W 

Instead of performing matrix transformations the stiffness, [K] , and mass, [M] , 
matrices are obtained using the La Grange formulation which states: 


and 
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where the potential energy, V, and the kinetic energy, T, are: 


V = 



(7) 


T 



. 2 , 
u dx 


( 8 ) 


After evaluating the integrals, the stiffness matrix produced by the new 
method is: 
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sym 


The mass matrix is: 



S 


v 2 
xm r 
2 



Since the first row corresponds to the displacement at x = 0, the boundary 
constraint, = 0, requires that the first row and column be deleted for 
calculation of the cantilever modes. 

The results of the modal solution are the frequencies and generalized dis- 
placements £ k = 1 , 2 .. . . The actual mode shapes may be obtained 

ox lei 

from equation (4). However, a mode acceleration method (U1MPROVE) available 
in NASTRAN and also used by Rubin will enhance the vectors. Transforming 
the matrix equations into equivalent integrals results in the equation: 
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where u(x) is obtained by multiplying the displacement u(x) in Equation (4) 
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by -co . This results in the mode shape: 
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where x = x/E is used 


lor simplicity. 


( 12 ) 


The exact solutions for the cantilever rod problem modal frequencies are: 


(2n - l)jf ft 
““ “ 21 *p 


ex 


The exact mode shapes are: 
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- sine S2R _ =M J 12£ 
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The calculated natural frequencies for the synthesized system produce an error 
ratio e , defined by the equation: 


c. 
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(15) 


The resulting errors in natural frequency are tabulated in Table 1. These 
errors match Rubins results for his method exactly. Note that, except for 
the single degree of freedom problem, the error in the last mode for any 
matrix size is nearly constant and that the convergence rate for a given 
order matrix is nearly uniform. 


An order of magnitude fit of the frequency errors is produced by the equation: 


e n. 
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where is Che frequency of Che lowest truncated mode shape. The equation 

is not accurate for the lower modes of the large order matrices due to compu- 
ter round-off. Single precision arithmetic produced numerical errors of the 

order 10 ^ * 


The RMS errors for the calculated eigenvectors are shown in Table 2 for the 
same order synthesized matrices. Both first and second methods for calcu- 
lating the vectors were used. The equations used for the vector errors are: 


E 


1 



<+i ' V dx 


(17) 


and 



Ct|, 2 " ^x )2dx 


(18) 


Both vectors were normalized to unit modal mass, 

Note that the improved displacement calculations (t|>„) produce much better re- 

« 5 

suits when the first order errors are between 10 and 10 

In other words, a good first approximation will produce a better improved so- 
lution. A poor first solution, such as the last mode in a set, will result 
in little improvement. A nearly exact first solution will not improve due 
to numerical round-off. 

The results of this test are nearly identical to Rubin's [Ref. 3] results 
for his method. The frequency errors fall exactly on the published curves. 
The displacement errors for the UAI method appear to be better than Rubin's 
results. However it is suspected that differences in numerical procedures 
produced these changes. Also the first order displacement results compare 
with the referenced results for the modified Bamford method used in Ref- 
erence 3. 
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TWO COMPONENT TRUSS PROBLEM 


This problem has nearly become a standard for the evaluation of modal syn- 
thesis methods. It has been used in References 4, 5, and others for com- 
parison between different formulations and procedures. A large quantity 
of data is available for validation of any new method. The basic problem, 
shown in Figure 2, consists of two truss substructures. Each nubstructure 
is reduced to its normal modes plus any additional shape functions used by 
a particular method. The trusses are combined at the three common grid 
points and the unconstrained modes of the combination are obtained. 

This problem was solved directly on the UAI modal synthesis system Implemented 
on L16 NASTRAN, Several different options and matrix sizes were tested. The 
parameters of the test cases are shown in Table 3. The matrix sizes were 
chosen to provide direct comparison with the results in Reference 4. 

The results were compared with a single-structure NASTRAN execution to obtain 
the percentage errors of the frequencies. These errors are shown in Tables 4 
and 5 along with results from HIntz [Ref. 4]. In all cases an excellent cor- 
relation was obtained between the NASTRAN results and the results of the 
equivalent formulations used by HIntz The only deviation occurred when the 
errors became too small to calculate when the NASTRAN printout truncated the 
difference in results. In each case the results are not shown where the 
eigenvector became unrecognizable and/or the natural frequencies changed in 
sequence. 

It Is important to note that the cases using free component: modes, with no 
inertia relief effects, produced very poor results. This is due to the fact 
that the free modes approximate half waves while the cantilever modes appro- 
ximate quarter waves. The shapes of the first modes of the combination 
apparently are difficult to approximate by a set of higher order shapes. 

The inertia relief shapes supply these smooth functions. Their contribution 
is most significant in the lower frequency modes. 

Also indicated in the tables, by dashed lines, are the lowest truncated fre- 
quencies for the component modes used in the analysis. The results indicate 
that in the inertia relief cases, this frequency is a good indication of the 
limit for valid results. Wien only normal and constraining modes are used, 
this frequency lias some significance, but does not Indicate possible errors 
due to poor approximation of the actual mode shapes. 

As a further check the problem was executed using 36 elastic degrees of free- 
dom (case 9). This case also correlated with the Reference 4 results, 
having 29 modes with a frequency er'ror of less than 5% Nearly all of the 
first 15 modes were calculated to values exact to the last digit of the print- 
out. However this case should not be considered ay a typical example since 
only 60 degrees of freedom existed in the original structure. The typical 
application of modal synthesis would result in a matrix size with a much 
smaller fraction of the original matrix size. 
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CONCLUSIONS 


As was postulated In Reference 1, the new modal synthesis method to be avail- 
able in NASTRAN is capable of simulating the results and accuracy of any of 
the current state-of-the-art modal synthesis methods. The differences in the 
results occur from selecting different types of component modes and types of 
solution vector recovery processing. Furthermore, it was observed from these 
tests that the frequencies of the truncated component normal modes are a 
significant indicator of the upper limit of valid combination modes. 

Although the new system provides accuracies equal to or better than any other 
advanced method, It also eliminates the restrictions that are Imposed by the 
other formulations. The UAI method does not require unconstrained modes 
required by the Rubin and McNeal [Rof. 6] formulations. The method conven- 
iently uses the actual boundary grid points as degrees of freedom (as in the 
Rubiti and MacNeal methods) as op >sed to the conventional Rayleigh-Ritz 
methods, in which an actual boundary displacement coordinate must be expressed 
as a combination of mode displacements. Furthermore, it allows any choice of 
mode shapes, including modes fixed at non-boundax"y points, partially free 
modes, and user supplied vectors. 

The results for both test problem ■ indicate that the inertia relief option 
is recommended for most cases. Trie number of calculations to obtain these 
shape functions is small relative to the modal calculations* A maximum 
of six extra degrees of freedom per component are added to the system. 

Results from the second test problem indicate that one should not replace 
modal coordinates with the inertia relief components since this will lower 
the effective frequency range. 

The use of the "improved displacement" options in the solution vector re- 
covery nrocess appears to be less dramatic in its effectiveness. This option 
will be most effective when the first order vectors ..re reasonably valid and 
accurate stress data are required. 
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TABLE 1. MODE FREQUENCY ERRORS, EW, VERSUS 
MATRIX ORDER - FREE ROD PROBLEM 
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TABLE 2. EIGENVECTOR RMS ERRORS VERSUS 
MATRIX ORDER - FREE ROD PROBLEM 


MODE 

NO. 


MATRIX ORDER 


1.31-2 3.18-5 3.06-6 6.16-7 1.81-7 1.35-7 2.23-7 

(2.69-2)* (4. 70-4) (9.76-5) (3.25-5) (6.97-6) (2.32-6) (6.52-7) 


9.81-3 1.03' 3 1.66-5 1.87-5 9.71-6 5.00-6 

(2.78-2) (3.79-3) (1.05-3) (2.01-4) (6.72-5) (1.25-5) 


2.15-2 3,36-3 2.16-4 5.44-5 5.50-5 

(4.30-2) (7.81-3) (1.11-3) (3.27-4) (6.00-5) 


3.19-2 

1.65-3 

2.77-4 

1.11-4 

(5.37-2) 

(4.10-3) 

(1.039-3) 

(1.73-4) 


9.19-3 1.12-3 2.04-4 

(1.51-2) (2.76-3) (4.05-4) 


4.79-2 

4.06-3 

5.46-4 

(6.79-2) 

(8.91-3) 

(8.34-4) 


1.56-2 

9.55-4 


(2.12-2 

(1.62-3) 


6.01-2 2.01-3 

(1.34-2) (3.06-2) 


3.96-3 

(5.96-3) 


9.12-3 

( 1 . 21 - 2 ) 


2.46-2 

(2.93-2) 


7.57-2 

( 8 . 20 - 2 ) 


*( ) without U IMPROVE 





















TABLE 3. .TEST CASE PARAMETERS 
FOR THE TWO COMPONENT TRUSS PROBLEMS 


CASE 

NUMBER OF 
MODAL COORDINATES 

NUMBER OF INERTIA 
RELIEF COORDINATES 

TYPE OF 

COMPONENT MODE 

TOTAL 

EuASTXC 

DOF 

1 

9 

0 

Free 

12 

D 

3 

6 

Free 

12 

■ 

17 

0 

Free 

20 

H 

11 

6 

Free 

20 

5 

9 

0 

Cantilever 

.12 

6 

3 

6 

Cantilever 

12 

7 

17 

0 

Cantilever 

20 

8 

11 

6 

Cantilever 

20 

5 

27 

6 

Cantilever 

36 
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TABLE 4. PERCENT FREQUENCY ERRORS WITH 12 ELASTIC DEGREES OF FREEDOM 

TWO COMPONENT TRUSS PROBLEM 


ELASTIC 

MODE 

NO. 

NASTRAN CASE 

REF. 

4 RESULTS 

1 

Free 

Modes 

2 

Free 

w/l.R. 

5 

Cant . 
Modes 

6 

Cant. 

v/l.R. 

Free 

w/I.R. 

Hurty 

Cant. 

w/I.R. 

J. 

24.27 

.006 

.011 

.00039 

.006 

.011 

.00067 

2 

3.28 

.021 

.013 

^194_ 

.019 

.013 

.187 

3 

10.41 

.737 

.031 

.737 

.074* 

.031 

.743 

4 

4,51 

.147 

.150 

2.93 

.150 

.155 

2.94 

5 

2.47 

1.82_ 

.197 

10.83 

1.68 

.190 

10.3 

6 

4.50 

6.45 

.184 

17.06 

6.55 

.184 

16.9 

7 

1.00 

16.02 

6.49 


16.8 

7.39 


8 

4.87 


6.44 





9 

0.75 







10 

11 


"N^indicates freq. of 
first truncated 
mode 





*Suspected typo error 
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TABLE 5. PERCENT FREQUENCY ERRORS WITH 20 ELASTIC DEGREES OF FREEDOM 

TWO COMPONENT TRUSS PROBLEM 



NASTRAN CASE 

REF. 

4 RESUTTS 

bLAo 1 1U 

MODE 

NO. 

3 

Free 

Modes 

4 

Free 

w/l.R. 

7 

Cant . 
Modes 

8 

Cant . 
w/I.R. 

Free 

w/I.R. 

Hurty 
(Cant . 
Modes) 

Cant . 
w/I.R. 

1 

8.92 

.00034 

.00043 

.00034 

.000017 

. 00074 

9xl0~ 9 








-6 

2 

1.21 

.00902 

.0017 

0 

.000061 

.0018 

3x10 

3 

7.67 

.0135 

.0098 

.0061 

.0138 

.0096 

.00584 

4 

1.08 

.00023 

.0096 

.00002 

.00024 

.0092 

.00002 

5 

6.00 

.00083 

.033 

0 

.00081 

.034 

.0014 

6 

0.85 

.0020 

.0098 

.00060 

.0020 

.0103 

.00054 

7 

0.61 

.080 

.947 

.268 

.083 

.941 

.264 

8 

1.58 

.0071 

.122 

.021 

.0068 

.117 

.018 

9 

.084 

.00098 

.59 

.54 

.00093 

.80 

.69 

10 

.030 

.0041 

.36 

.40 

.0045 

.20 

.25 

11 

.90 

.021 

.33 

_ . 98_ 

.022 

.30 

1.03 

12 

3.30 

.428 

.49 

12.3 

.134 

.28 

11.1 

13 

4.01 

5 . 35_ 

.16 


5.33 

.14 


14 

.244 

7.87 

.77 


7.15 

.72 


15 

1.10 


2.37 



2.63 


16 

-.59 


12.15 



11 . 4 


17 

6 . 69_ 


mm> 







< sT indicates freq. of 






^first truncated mode 
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ORIGINAL PAGE IS 
OP POOR QUALITY! 



(a) Basic Definition 




(b) Component Modes (Free) tj> 


= cos 


mrx 
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(c) Solution Modes (f) 


sin 


(2n - 1)ttx 
2 £ 


FIGURE I. EXACT ROD PROBLEM 




o 


• - Boundary Points 



Figure 2. Nine-cell truss basic substructures. 





